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ABSTRACT 

We demonstrate that the corotation resonance affects only some non- 
axisymmetric g-mode oscillations of thin accretion disks, since it is located within 
their capture zones. Using a more general (weaker radial WKB approximation) 
formulation of the governing equations, such g-modes, treated as perfect fluid 
perturbations, are shown to formally diverge at the position of the corotation 
resonance. A small amount of viscosity adds a small imaginary part to the eigen- 
frequency which has been shown to induce a secular instability (mode growth) 
if it acts hydrodynamically. The g-mode corotation resonance divergence disap- 
pears, but the mode magnitude can remain largest at the place of the corotation 
resonance. For the known g-modes with moderate values of the radial mode 
number and axial mode number (and any vertical mode number), the corotation 
resonance lies well outside their trapping region (and inside the innermost stable 
circular orbit), so the observationally relevant modes are unaffected by the reso- 
nance. The axisymmetric g-mode has been seen by Reynolds Sz Miller in a recent 
inviscid hydrodynamic accretion disk global numerical simulation. We also point 
out that the g-mode eigenfrequencies approximately obey the harmonic relation 
a oc m for axial mode numbers \m\ > 1. 

Subject headings: accretion, accretion disks — black hole physics — gravitation 
— relativity 
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Introduction 



In principle, all adiabatic perturbations of equilibrium models of accretion disks can 
be analyzed in terms of global normal modes. The pioneering studies of Shoji Kato and 
his group and the more recent work of our group have focused on accretion disks around 
black holes, so that no complications fr om boundary layers are involved. For a recent re- 
view of 'relativistic diskoseismology', see iKatol (120011 ). A short summary of observationally 



relevant results from our recent analyses of the low-lying spectrum, which consists of g 



modes (jPerez. Silbergleit. Wagoner fc Lehr 



( iSilbergleit. Wagoner fc Ortega- Rodriguez 



1997) (hereafter referred to as RD.I), c-modes 



2001) (hereafter refer red to as RD.II), funda- 



mental p-modes ( Ortega- Rodriguez. Silbergleit fc Wagoner 2002^ fhereafter referred to as 
RD.III), and other p-modes (lOrtega- Rodriguez. Silbergleit fc Wagoner! 120071) fhe reafter re- 
ferred to as RD.IV), is given by lWagoner. Silbergleit. fc Ortega- Rodriguez! (120011 ) . 



Local analyses (restricted radial interval) have also played an im portant role in our 
understanding of these perturbations ( Kato. Fukue fc Mineshigel 1998). In deed, one of the 
motivations for this paper was the search by lLi. Goodman fc Narayanl (120031 ) for a dynamical 
instability via such an analysis of traveling waves impinging on the corotation (and Lind- 
blad) resonance. They argued that no modes were likely to be dynamically unstable, since 
amplification at the corotation resonance would not occur in realistic thin disks. 

In section 2 we summarize the foundations (assumptions and equations) of our approach. 
In section 3 we examine the corotation resonance location and show that only certain non- 
axisymmetric g-modes may be affected by the resonance. In section 4 we investigate the 
behavior of the vertical and radial eigenfunctions near the corotation resonance, exhibiting 
the local and global divergence, and thus prove that the range of the eigenfrequencies of non- 
axisymmetric g-modes is reduced to its upper part (specified by the rotational frequency, 
radial epicyclic frequency, and the azimuthal wave number of the mode). In section 5 we 
discuss effects of introducing viscosity and buoyancy. We also comment on the possible 
observational relevance of the spectrum of the g-modes. 



2. Basic Assumptions and Equations 



We take c = 1, and express all distances in units of GM/c 2 and all frequencies in units 
of c 3 /GM (where M is the mass of the black hole) unless otherwise indicated. We employ 
the Kerr metric to study a thin accretion disk . The equilibrium disk is taken to be describe d 
by the standard relativistic thin disk model (INovikov fc Thorndll973l ; iPage fc Thorndll974l ). 
The velocity components v r = v z = 0, and the disk semi-thickness h(r) ~ c s (r,0)/Q <C r, 
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where c s (r, z) is the speed of sound. The key frequencies, associated with free-particle orbits, 
are 

Q(r) = (r 3/2 + a)' 1 , 
Sl ± (r) = fi(r) (1 - 4a/r 3/2 + 3a 2 /r 2 ) 1/2 , 
k(t) = fi(r) (1 -6/r + 8a/r 3 / 2 -3a 2 /r 2 ) V2 (1) 

the rotational, vertical epicyclic, and radial epicyclic frequencies, respectively. The angular 
momentum parameter a = cJ/GM 2 is less than unity in absolute value. The inner edge of 
the disk is at approximately the radius of the last stable free-particle circular orbit r = Tj(a), 
where the epicyclic frequency «(rj) = 0. So all the relations we use are for r > r,, where 
n(r) > 0. The outer disk radius is denoted by r Q . 



We apply the general relativistic formalism that llpser fc Lindbloml (119921 ) developed for 
perturbations of purely rotating perfect fluids, although the effects of viscosity are included 
in the equilibrium model (and some of our notation differs from theirs, see below). The 
pressure p is much less than the mass-energy density p. We neglect the self-gravity of the 
disk, which is a good approxi mation since the rat i o of d isk to black hole mass is usually very 
small (see RD.IV, section 7). llpser fc Lindbloml (119921 ) then show that one can express the 
Eulerian perturbations of all physical quantities through a single function 



8V 



5p 

p(3u 



(2) 



Due to the stationary and axisymmetric equilibrium, the angular and time dependences are 
factor ed out as 5V = V(r, z) exp[i(m(j)+at)]) , and the master equation (39) of llpser &: Lindblom 
(119921 ) for the function V [see also RD.I, equation (2.21)] assumes the form 



ld_ 

r dr 



r 2 T 3 p 



UJ — K 



dV 

2 I Q r 



+ 



d_ 

dz 



rTp 



uj- 



lu 2 -N 2 dz 



dV 



+ rTf3uj^V = 



The corotation frequency to is related to the eigenfrequency a by 

u)(r, a) = a + mQ(r) . 



(3) 



(4) 



The buoyancy frequency is dominated by its vertical component N z (RD.I). The Kerr metric 
component g rr = T 2 (r) and four- velocity component u° = dt/dr = (3{r) are both of order 
unit y (in Boyer-Lindquist co ordinates); their expressions can be found in RD.I. Compared 
with llpser &: Lindbloml (119921 ) . the definitions of a and to are switched, while their 7 becomes 
our p. The function $(r, z, u,m) depends on various properties of the unperturbed disk, as 
well as the mode eigenfrequency, a, and its azimuthal wave number, m. 
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To simplify the analysis, we consider barotropic disks [p = p(p)], so the buoyancy fre- 
quency N z = 0. Since this frequency should be lower than other characteristic frequencies, 
and in general vanishes on the mid-plane of symmetry z — 0, this should be a good approx- 
imation except possibly in the neighborhood of the corotation resonance [where ou(r) = 0]. 
Hydrostatic equilibrium then provides the following unperturbed density and pressure dis- 
tributions (r > 1 is the adiabatic index): 

p(r,y) = p (r)(l-y 2 y , p(r, y) = p (r)(l - y 2 ) 9+1 , g = 1/(T - 1) > ; (5) 

the disk surfaces are at y = ±1, with the new coordinate y related to the vertical coordinate 
z through the characteristic disk semi-thickness h(r): 

- J— l T ~ 1 hi \ - 1 / Po( r ) ^ 

V h(r) V 2r ' n[r) /3(r)^ ± (r) y Po(r) ' W 

The speed of sound c a (r,y) is specified by 

c 2 = r J 9/ P = r(^) 2 (i- z/ 2 ). (7) 



With these simplifications, for the function $ which appears at the end of the master 
equation ([3]) we obtain: 



$ (r,z) = ^- + 



5 



pT 3 r 2 9 2pTrfi 2 (m + 



Tr 2 dr \_[3 2 {uj 2 — k 2 ) dr f3 3 (u 2 — k 2 ) 

The quantity Q z is a component of the angular velocity four-vector, while is another 
component of the fluid four-velocity. Like all other quantities in the expression (JBJ) except p 
and c s , they are functions of r alone. 

Because of the structure of the density and speed of sound as functions of the coordi- 
nates, the equation ([3]) does not allow for an exact separation of variables in either of their 
two pairs, {r, z) or {r, y}. Therefore in the past we adopted the (usually realistic) assump- 
tion of strong variation of the modes in the radial direction (characteristic radial wavelength 
A r <r), and used the asymptotic separation of variables based on it. In particular, we look 
for a separated solution to the master equation flS} of the form 

V = V r (r)V z (z,r) (9) 

Unlike the previous analyses, here we do not neglect the derivatives of functions of u(r) [in 
addition to V r and (uj 2 — k 2 )^ 1 ], since it varies much more strongly than V z and all other 
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quantities near the corotation resonance (for m/fl). However, we take all other functions 
out of the radial derivatives. Then, with the separation (Q also used, the master equation 
assumes the form 



X 2 d 

V r dr 



UJ 2 — K 2 



dV r 
dr 



+ G{r) 



1 



d 



2\g 



OK 
dz 



to: 



We have introduced ujJr) 



{l-ytyV z dz 
u>/£l±, and the function G(r) is defined by 

c. 



I7i 2 (l -y 2 ) 



G(r) = G(r, a, a, m) = 



P 



(10) 



(11) 



T 2 u 



d_ 

dr 



K 



did 

dr 



2Q z u d f m + l3u ( j ) u)\ (m + jSu^cu) 2 ^ 2 
j 2 -k 2 ) ~ (3 2 r 2 T 2 {u 2 - k 2 ) 



f3 2 r dr 

The left-hand side of equation ( flOl) is a rapidly varying function of r, while the right- 
hand side is a rapidly varying function of z. Therefore within this (weak) radial WKB 
approximation both are equal to a slowly varying separation 'constant' S(r), which (in 
keeping with our previous convention) we denote instead by 

S(r) = Wr)-uj*}/(Th 2 ). (12) 

The vertical eigenvalue ^ is thus the redefined separation 'constant'. Employing everywhere 
the new vertical coordinate y defined by equation ([6]), this radial WKB approximation then 
does indeed produce separated equations for V y = V z (y,r) (a slowly varying function of r) 
and V r (r): 



d 2 V 

d_ IV u 2 

dr \ to 2 — 



2gy 



dy 
dV r 
dr 



+ 2g [uly 2 + ^{l-y 2 )]V y 







+ G 



Th 2 



~£2 







(13) 



(14) 



Away from the corotation resonance, G ~ 1/r 2 + m 2 /r 2 , which is seen to be of order ihjr) 2 
smaller than its competing term u 2 /(Th 2 ) = [j3u/c s (r, 0)] 2 in equations ffTUl) and f[T4"|) . That 
is why G(r) has been neglected in RD.I-IV. 

Together with the proper boundary conditions, discussed in detail in the referenced 
papers, these two equations specify (slowly varying) vertical eigenvalues \1/ = \l/(cr, r) and 
eigenfrequencies a, as well as the corresponding vertical (V y ) and radial (V r ) eigenfunctions 
for modes of all types. Lagrangian displacements are related to V by the expressions 



uT 2 



dV 



(5{ui 2 — k 2 ) dr 



1 dV 

(3uj dz 



/3 2 [l-(2/r)(l-a^)] 



ltd 



dfl ru z 
dr j3 2 (r 2 — 2r + a 2 ) 



(15) 



[See equations (2.25) and (2.26) of RD.L] Here u z is a component of the vorticity four-vector 
jlpser & Lindblomlll992h . 
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3. The Corotation Resonance and Its Implications for Diskoseismology Modes 

3.1. Frequency range and location of the corotation resonance 

As always, in view of the {cr — > —a, m — > — m} symmetry, we can restrict our consid- 
eration to m > 0. For axially symmetric modes m = and u = a ^ 0. However, for m > 
there might exist a point, r cor = r cor (a, m, a), where the corotation frequency ui goes to zero, 

u{r cor ) = a + mfi(r„) = . (16) 

This specifies the location of the corotation resonance for a given eigenfrequency. Since Q(r) 
is a decreasing function of the radius, the unique solution r cor = r cor (a,a,m) to equation 
( I16p exists within the disk, < r cor < r Q , for all eigenfrequencies in the range 

mQ(r ) < —a < mO(rj) . (17) 

It is important to compare r cor with the Lindblad resonances at the radii r± = r±(a, a, m) 
(r_ < r+) defined as the roots of u 2 (r) — k 2 {t) = 0. As pointed out in RD.I (note especially 
fig. 3), for m > they exist in a wider range of the eigenfrequencies than (fTTl) . namely, when 

mQ(r ) < —a < max [mQ(r) + «(r)] . (18) 

Ti<r<r 

This is also the maximum possible frequency range of g-modes, according to their definition 
specifying their capture zone as the interval between the Lindblad resonances where /t 2 (r) — 
uj 2 (r) > 0. 

So, whenever the corotation resonance exists, the Lindblad resonances are also present 
(but not necessarily vice versa). Moreover, since n(r) is positive inside the disk, we can write 

mQ(r_) = —a + /t(r_) > —o = mQ(r cor ) > —a — k{t + ) = mQ(r + ) 

Note that the inequalities here are based on the expressions involving a . By looking at their 
counterparts with the angular velocity Q(r) and invoking again its monotonic decrease, we 
conclude that 

r_ < r cor < r+ . (19) 
Hence, whenever the corotation resonance occurs, it lies between the Lindblad resonances. 

3.2. Corotation resonance and the diskoseismology modes 

The result ([19]) shows that the corotation resonance does not significantly affect both 
c-modes, captured in a region < r < r c < r_ of the inner disk region (see RD.II), and 
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inner and outer p-modes, residing respectively in < r < r_ and r + < r < r Q (see RD.III, 
IV). The only disk oscillations which could be strongly affected by it are thus the (non- 
axisymmetric) g-modes, since the region r_ < r < r + is their capture domain. So, the first 
significant result of the paper allows us to discuss only g-modes in the sequel. It should 
also be noted that g-modes are the most robustly determined, since their capture zone does 
not include either of the disk boundaries, where the physical conditions and validity of our 
assumptions are more uncertain. 

Furthermore, for the effect of the corotation resonance to show up, the resonance itself 
must be present; this means that only those g-modes can be influenced by it whose eignefre- 
quencies are in the corotation resonance range (fTTj) . All the g-mode eigenfrequencies found 
so far (see RD.I, Tables 1 — 3) belong, in fact, to the more restricted range 

mQ(ri) < —o < max [mQ(r) + n(r)} , (20) 

ri<r<r 

i.e., to the upper part of the maximum g-mode range (1181) . The next section shows that this 
is not unexpected, since the g-modes with the eignefrequencies in the range (fTTj) are simply 
absent. 



4. Corotation Resonance Singularity and 
the Frequency Range of Non— Axisymmetric g— Modes 

4.1. Possible g— mode divergence at the corotation resonance 

Let us find out what can happen to those m > g-modes which have the corotation 
resonance in their capture domain, that is, eigenfrequencies in the interval (fTTj) - The behavior 
of the eigenfunction, V, and its radial derivative, dV/dr, at the corotation resonance needs 
to be examined, because one or both of them may become singular there. This would imply, 
in the first place, the singularity of the Lagrangian displacements by the formulas (fTBT) . which 
is hardly acceptable for small perturbations. However, as shown in section 4.3, this is exactly 
what happens. Unfortunately, an integral divergence turns out also to be involved. 

To see this, one multiplies the master equation (j3J) by rV(r, z), integrates the result over 
the disk, and then integrates by parts taking into account that the unperturbed density (p) 
vanishes at the disk boundary, to obtain: 

J r 2 Tf3u§ V 2 drdz = J r 2 Tp 

disk disk 



dV 

dz 



T 2 co 2 



+ 



fdV 

\ dr 



drdz 



(21) 
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The left hand side of this identity must be finite for any reasonable perturbation, since it 
represents the relevant norm of it. The two sources of divergences on the right are: (a) the 
corotation resonance, and (b) the Lindblad resonances. 



For the source (b), the combination 



K 



2\-l 



(dV/dr) has been shown to be finite 



at r± in RD.I and RD.III, and this conclusion holds with some new singular terms in the 
coefficient G(r): by its definition (TTTj) . G(r) = O ((u 2 — k 2 )~ 2 ) when r — > r±. Therefore 
the second term on the right of the identity (1211) is finite (in fact, vanishes) at the Lindblad 



resonances. Moreover, V and {uj 2 



K 



L V are both finite there as well, so that there is no 



Lindblad resonance divergence in either the first term on the right of equation (121]) . or in 
the term on the left stemming from the singularity in G(r). 

However, it turns out that at the corotation resonance, the l.h.s. remains finite but 
both terms on the r.h.s. of the identity f[2~Tl) generically diverge. To show this, we start with 
the first term, and transform it using the WKB separation, redefinition (jSj) of the vertical 
variable, and expression dSJ) for the density: 



disk 



r 2 Tp (^—) drdz 
\dz J 



r 2 Tp 



k 



2T 



V 2 dr 



2\9 



-1 



dVy 

dy 



dy 



V 2 L, dr . 



(22) 



Here k > is the proper average value of the positive function r 2 Tp h 1 a/(L — 1)/2T of the 
radius, and we denoted 



I v (r)= I ^-y 2 ) 9 [Vy{r,y) 



dy 



(23) 



(Prime will now denote the derivative with respect to y.) Clearly, I y {r) > for any r; 
if I y (r*) = for some r = r*, then equation ( |23l) implies V y (r*,y) = 0, so the vertical 
eigenfuntion V y (r*,y) = const. However, a constant never satisfies the vertical equation fTTSl 
except for the case r = r cor , ^(r cor ) = 0. The last condition is never true for g-modes; 
the demonstration requires nevertheless a slightly more sophisticated argument (following 
immediately) . 



4.2. Vertical eigenvalue problem near the corotation resonance 



We rewrite the vertical equation fflBl in the self-adjoint form, 



+ 2g 



(i-y 2 y- 1 v y = o 



(24) 
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multiply by V y and integrate over the interval (—1, 1). Using integration by parts and the 
fact that V y (±l) must be finite, we thus obtain 

Iy(r) = j\(l-y 2 ) 9 (vtfdy = 2g j\ 2 + $ j\l-y 2 ) 9 V 2 dy 

At the corotation resonance this equality becomes 

I y (r cor ) = 2gV(r cor ) J (l - y 2 ) 9 [V y (r cor ,y)] 2 dy ± , (25) 

unless \l/(r cor ) = 0. However, if this is the case, we can expand, in the vicinity of r cor , the 
eigenvalue and eigenfunction in the small parameter uj 2 = uj 2 /VL 2 ^ as 

V = V 1 ul + ..., V y = V^ + V {1) uol + ... , 

obtaining thus from the equation (|24l) : 

CV<®=0, CV^ = -2g [1 - (1 - y 2 )] (l-^f 1 ^...; = [(l-y 2 )V 

The only zero-order solution bounded together with its derivative at y = ±1 is, of course, 
V^(y) = 1, and the solvability criterion of the problem for the first correction, 

= £ V®CVW dy = -2g £ [l - ^ (l - y 2 )} (l - y 2 ) 9 ' 1 dy , 

provides = — (r — l)/2 by means of an easy calculation. Hence 

^(r)/cu 2 (r) = _(r-l)/2+0 (cu 2 (r)) = -(r~l)/2+0 ((r - r cor ) 2 ) < 0, r near r cor , (26) 

which, by definition, corresponds to some p-mode, and not a g-mode (characterized by 
/uj 2 > 1). Note that the last term in the equalities (126]) holds due to the fact that u, and 
hence u;*, is to lowest order linear near r = r cor : 

W (r) = m^^(r - r cor ) + O ((r - r cor ) 2 ) . (27) 



So \l/(r cor .) 7^ for any g-mode, and therefore I y (r cor ) 0. Thus, from equation ([22]) . 
the question of whether the first term on the right of ( 12T1) diverges or not at the corotation 
resonance reduces to the same question about the norm of the radial eigenfunction (whose 
square we denote I r ): 

/ bv\ 2 r° 

r 2 Tp — drdzocI r , I r = / V 2 dr . (28) 



dz J 



disk 

To obtain the answer, it remains only to investigate the behavior of V r (r) near r = r cor . 



lr The same derivation of this result was carried out in RD.III, section 3, for m = 0, cj* = er/i7j_; which 
specifics play, in fact, no role in it. 
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4.3. Behavior of the radial eigenfunction near the corotation resonance 



Using the Taylor expansion fl27|) . we see that the radial equation (fT4"l) near r cor can be 
written as 



d 2 V r 
dr 2 



+ 



r — r r 



where 



+ 0(1) 



Ql 



dVr 

dr 



+ 



[r — r r 



+ o(( 



r — r r 



o- 1 ) 



V r = 



/3Vni* fdn\~ 2 



— 



(29) 



(30) 



m 2 c 2 (r, 0) \dr J 

and all the higher order terms are integer powe rs of (r — r cor ). According to the analytical 
theory of second order ODEs [e. g., lOlverl (119821 )]. the general solution to equation fl29l near 



r cor has the form: 

V r (r) = C + (r - r cor r+(l + ...) + C_(r - r cor )^(l + ...). (31) 
Here C± are some constants, and i>± are the roots of the characteristic equation 



u 2 + u + Q 2 c 







(32) 



Note that for an isothermal thin disk model, iLi. Goodman Narayanl (120031 ) have found a 
different behavior of traveling waves. 



The results f l3~Tj) and (132^) show that in the vicinity of the corotation resonance 

V r (r) = 0((r- w)- 1 / 2 ) , dV r (r)/dr = O ((r - r cor )- 3/2 ) , r r cor . (33) 

The singularity is, in fact, even stronger if Q 2 < 1/4, but this is not important here. It is 
also rather an exceptional shown in RD.I, for any axial mode number rn > 1 and 

radial mode number n > 0, the vertical eigenvalue ^ = ^-(r) grows indefinitely with the 
vertical mode number j for all relevant radii, including r = r cor . So, by the formula ( J30l) . 
Q 2 > 1/4 for an infinite set of vertical eigenvalues, with the opposite inequality valid perhaps 
just for a few of them, if at all. Finally, the only case of a singularity weaker than the one 
in equations (1331) is Ql < 1/4 and C_ = in expression (j3~Ij) . However, this requires fine 
tuning of parameters. Therefore the estimates (|33|) always hold, as well as 



V 2 

< ... 



O ((r - rw)" 1 ) , uj 2 (dV r /dr) 2 = O ((r - rw)" 1 ) , r -> r c 



(34) 



implied by them. 
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4.4. Local and global corotation resonance divergences. 
Range of non— axisymmetric g— modes 

Formulas (133]) combined with the expressions (1T51) show that all the three components 
of the Lagrangian displacement are singular at the corotation resonance: 

f -(r-rw)- 1 / 2 , T, e ~ {r - r cm )-" 2 . 

However, note from equation (J2J) that the perturbation <5p/p oc (r — rw) 1 / 2 vanishes there. 

In addition to this local singularity of the physical quantities, a global one also exists. 
According to the first of the equalities (1341) . the integral I r defined in (1251) always diverges, 
hence so does the first term in the expression (l2Tj) for the norm of the eigenfunction. The 
second term there, in fact, also diverges via the second of the equalities ( 134)) . Such a com- 
bination of local and global perturbation singularities is definitely unacceptable. (The norm 
of an eigenfunction belonging to the discrete spectrum must be finite.) 

We have thus proved that non-axisymmetric g-modes with eigenfrequencies in the coro- 
tation resonance range (1T71) cannot exist within the framework of inviscid perturbations, so 
their actual frequency range (12"U1) . 

mVL{ri) < —a < max [mQ(r) + ft(r)] = ft(r m ) + mf2(r m ), m — 1, 2, . . . , (35) 

Ti<r<r 

is the upper part of the maximum possible range (1T5|) . for which the corotation resonance is 
absent within the disk. In contrast with this, the range of the axially symmetric g-modes, 

n{r ) < \a\ < max n(r) = K (a) (36) 

ri<r<r 

includes low frequencies as well. 



5. Discussion and Conclusions 

5.1. Non— axisymmetric g— modes with moderate radial and azimuthal wave 

numbers 

The established corotation resonance divergence could cast a shadow of doubt on the 
known results for non-axisymmetric g-modes, in particular, on those found in RD.I, espe 



daily since the WKB technique used there [see also lOrtega- Rodriguez. Silbergleit fc Wagoner 



(120071 )] to calculate eigenfrequencies is not fully sensitive to the presence of the corotation 



resonance. However, as pointed out in section 3.2, all the found eigenfrequencies belong to 
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the proper range fl35l) . \a mn j\ > mfifo) (m, rt, and j are the azimuthal, radial and vertical 
mode numbers). That means that the corotation resonance lies inside the inner radius of 
the disk (where the gas is spiraling into the black hole). We now indicate why the g-modes 
(with moderate values of m and n) are the most robust modes, and therefore astrophysically 
the most relevant. 

For any azimuthal mode number m, we denote the largest eigenfequency \a mn j\ by 
a m = n(r m ) + mQ(r m ) . In the limit j — > oo, the trapping zone Ar = r + — r_ — > 0, with 
r + > r rn > r_. We now consider moderate values of m and n, but any value of the vertical 
mode number j. Then from equation (5.3) and Tables 1-3 of RD.I, it is seen that these 
g-modes occupy a very small frequency range below their maximum: 

Thus, as also seen from Fig. 3 and Fig. 5 of RD.I, it follows that 

\<?mnj\ ~ mfi(ri) (37) 

for moderate values of n and m > 1, with the approximate equality becoming quickly more 
and more exact as m grows. Therefore for a given black hole, the largest frequency splitting 
should be due to the azimuthal mode number m = 0, 1, 2, . . ., with the eigenfrequencies 
| cr | oc m for m > 1. We also note that with increasing values of m, the mode location 
r m — y r i an d its extent Ar — > 0. Then the mode leakage into the flow onto the black hole 
and the uncertainties in the physical conditions near become more important while the 
fractional modulation of the luminosity decreases. Therefore the low m modes should be the 
most robust and observable. 



5.2. QPO features: 3/2 and other integer frequency ratios 



In this connection, one should note the observational claims that some of the quasi-stable 
high frequency QPOs in black hole binary X-ray sources have frequency ratios close to 3/2 
( IMcClintock fc Remillardl 120041 ; iRemillard &: McClintockl 120061 ). This can be, in principle, 
explained by excitation of two (groups of) g-modes with m = 2 and m = 3, as suggested in 
RD.IV. 

Any such explanation for any QPO frequency implies a relation between the mass and 
angular momentum of the black hole in the corresponding X-ray binary. Indeed, using the 
general relation for the dimensional frequency / = 3.23 x 10 A (M & /M)\a\ Hz, from formula 
f l37|) one obtains 



M/M Q = 3.23 x 10 4 mfi(ri)// m (Hz) , m > 1 . 



(38) 
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Since 7*j = rj(a), fi( r «) = (r/ + a) - = jF(a) is a universal function of black hole angular 
momentum only. It is a monotonically increasing function, with 1) = 0.038 (extreme 
counter-rotation), .F(O) = 0.068 (non-rotating black hole), and ^(l) = 0.5 (extreme rota- 
tion). 



Let us compare the prediction of equation (I3"8l) with the data (IMcClintock fc Remillard 



2004 ) for the three binary black holes with 3/2 or 3/2/1 frequency ratios and measured mass. 
Using the lower limit T{a) > 0.038, we obtain M > 8.2M for GRO J1655-40 (/ = 300,450 
Hz), observed to have M = 6.3±0.3M o . Similarly for XTE J1550-564 (/ = 92, 184, 276 Hz), 
we obtain M > 13.3M , compared to the observed M = 9.6 ± 1.2M Q . There are two pairs 
of 3/2 QPOs observed in GRS 1915+105 (/ = 41,67; 113, 168 Hz), whose black hole has a 
measured mass of M = 14 ± 4M Q . For the higher frequency pair, we obtain M > 22M Q 
while for the lower frequency pair we obtain M > 58M . If we use no rotation rather than 
counter rotation to provide the lo wer mass limits, they would inc rease by a factor 1.8. Similar 
conclusions have been reached by lTassev &: Bertschingerl (120071 ) . 



It is perhaps not surprising that these QPOs cannot be explained as m > 1 g-modes, 
since they occupy a much smaller r egion of the disk nearer to its inner edg e (where our 
accretion disk model is most suspect) (IPerez. Silbergleit. Wagoner fc Lehrlll997l ). In addition, 
their dependence reduces the observed modulation. 

One would thus wonder why the m = g- mode was not seen. It provides a relation 
MJM Q = F (a) I fo similar to equation (RD.I) (jWagoner. Silbergleit. fc Ortega- Rodriguez 
20011 ). Indeed, all but two of the QPOs in these sources could be fundamental g-modes 
(for some value of a < 1), but of course only one in each source. The value of a required for 
some of these is close to that estimated from the spec t rosco pic method [temperature and lu- 
minosity determine rj(a)]. For instance, IShafee et al.l (120061 ) obtain a = 0.65 — 0.75 for GRO 
J1655-40, whereas ident ifying its 300 Hz QPO as an m = g-mode requires a = 0.9 — 1.0. 
McClintock et al.l (l2006h obtain a > 0.98 for GRS 1915+105, whereas identifying its 168 Hz 
QPO as an m = g-mode requires a = 0.8 — 1.0. 



From their MHD simulations, iTagger fc Varnierd (120061 ) claim that the m = 2,3,... 
g-modes can grew to dominance over those of m — 0, 1 via the Rossby wave instability. 
However, this required a large concentration of magnetic field between the black hole and 
the accretion disk. I n addition, their simula t ion ne glected the vertical structure of the disk. 
We should note that lArras. Blaes fc Turner! (120061 ) found no g-modes in their shearing-box 
MHD simulations of a limited radial region of an accretion disk. There were some indications 
of the generation of p-modes within the MRI-induced turbulence, however. 



Very recently, 



Reynolds fc Miller 



(120071 ) reported results of ideal hydrodynamic (2D and 
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3D) an d MHD (3D) global numerical simulations of accretion disks. As in lTagger fc Varniere 
( 120061 ). the nonrotating black hole was represented by a modified Newtonian gravitational 
potential. The evolution was typically followed for about 10 2 orbital periods of the inner disk. 
From power spectra at many radii, the m = g-mode was seen in the hydro simulations 
at the predicted frequency and radial extent. It was not seen in the MHD simulations. 
However, because of the induced MRI turbulence, it would not be expected to be seen if 
it was at the same amplitude as in the hydro simulations. Because of the limited range of 
(with periodic boundary conditions) and frequency in the 3D simulations, the higher m 
modes could not have been seen. 



5.3. Perturbative effects of buoyancy and viscosity 



We need to say a few words about a non-zero buoyancy manifested by a Brunt- Vaisala 
frequency N z = N z (r,z) > 0. This frequency is involved in all the equations only via the 
expression {oj 2 — N 2 ), so the corotation resonance equation (FTO]) becomes 



uj 2 - n: 



o 



One might expect that the buoyancy would be small in re alistic accretion disks, s ince the 
magneto- rotational instability produces strong turbulence (IHawley &: Krolikl l200ll ). which 
should locally homogenize the specific entropy. If so, N z should be treated as a perturbation 
in the above equation. Because of that, it leads only to a small and generally z-dependent 
change in the corotation resonance position. In fact, the corotation resonance point splits 
into two nearby ones, 



/yi — ^— ryt 

' cor ~ c 



± 



iV (r 



m 



dVt(r c 



dr 



Our arguments and results regarding the mode divergence apply to both of them. 

In the presence of viscosity which acts hydr odynamically (via the a model) and per- 
turbatively, lOrtega-Rodrfguez fc Wagoner! (120001 ) found that for most (including these g-) 



modes, a viscous instability is induced. Such accretion disks are thus secularly unstable. 
An effective viscosity (in particular, generated by the magneto-rotational instability) should 
be present in these (thin) accretion disks, but it is not known in what ways it acts like a 
hydrodynamic viscosity. To lowest order in its magnitude, the viscosity does not change the 
values of the eigenfrequencies, but introduces a small i maginary part in them. As usual, this 
imaginary part removes the divergences found above (jNowak fc Wagoner! Il992l ). since the 
corotation resonance equality (Tl6|) no longer holds at any radius (a well-known effect of a 
complex pole near the real axis). However, the actual oscillation amplitude will usually still 
be largest at the corotation resonance radius as long as the mode amplitude is linear. 
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One is naturally tempted to contemplate the outcome of the mode growth. Will enough 
nonlinearity be induced to lead to significant mode- mode coupling and relate d effec ts? One 
would like to extend the local resonance analysis of lAbramowicz fc Kluzniakl (120011 ) to this 
problem. However, one must also investigate how much the spatial extent of the modes aver- 
ages over the strong MRI-induced turbulence. Could it be sufficient to provide a stochastic 
driving force (producing mode excitation) within a time-independent 'unperturbed' state? 



This work was supported by NASA grant NAS 8-39225 to Gravity Probe B. 
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